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Abstract. Fractional equations appear in the description of the dynamics of various 
physical systems. For Lagrangian systems, the embedding theory developped by Cresson 
['Fractional embedding of differential operators and Lagrangian systems", J. Math. Phys. 
48, 033504 (2007)] provides a univocal way to obtain such equations, stemming from a 
least action principle. However, no matter how equations are obtained, the dimension 
of the fractional derivative differs from the classical one and may induce problems of 
temporal homogeneity in fractional objects. In this paper, we show that it is necessary 
to introduce an extrinsic constant of time. Then, we use it to construct two equivalent 
fractional embeddings which retain homogeneity. The notion of fractional constant is also 
discussed through this formalism. Finally, an illustration is given with natural Lagrangian 
systems, and the case of the harmonic oscillator is entirely treated. 



1. Introduction 

For about twenty years, the fractional calculus has known significant development, ac- 
cording to its successes in various domains, such as chaotic dynamics, viscoelasticity, acous- 
tics, electricity or polymer chemistry [H El 13, El [lO] . In this context, fractional equations 
are used for modeling physical phenomena. No more considered as formal objects, but as 
physically significant, one may expect them to be homogeneous, a fundamental principle 
of the physical equations. In other words, the objects of a same equation must have the 
same dimension, expressible in terms of fundamental units, such as length L, mass M or 
time T. For example, the dimension of a force F is [F] = M L'^T~'^. 

However, in many fractional equations we can find in litterature, homogeneity is not 
respected. Let firm up details of our topic: from a strictly mathematical point of view, this 
notion of homogeneity is meaningless. The problem arises only if the equations considered 
are susceptible to model a physical system. But in that case, this question is fundamental. 
However, the problem is often hidden in modelisations because the constants are seen as 
parameters to fit. 

How to obtain fractional equations in time ? The usual approach consists in subsistuting 
in the classical equation considered the derivative d/dt (homogeneous to T~^) by a new 
one, d^ /dt^ (homogeneous to T~"), with < a < 1. Clearly, this transformation does not 
conserve the temporal homogeneity of the equation. This remark would apply to fractional 
derivatives in space, but in this article, only the temporal ones will be considered. 

The problem of the passage from the classical case to the fractional one has already been 
discussed in [Ij. Shifting the derivative leads to fractional equations, but, from a physical 
point of view, is it the right procedure ? Do the fractional equations obtained remain 
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physically relevant ? Answers are given in the case of Lagrangian systems. Extending the 
least action principle, Cresson [4] provides an univocal method, the fractional embedding, 
which leads to fractional equations. The respect of this other fundamental principle makes 
the equations susceptible to conserve a physical meaning. 

This article is in keeping with this philosphy. In section [21 it is shown that an extrinsic 
characteristic time has to be introduced. After having presented the notion of embedding 
in section [HI we show in section [H that this time constant leads to two equivalent methods of 
fractional embedding, which additionally conserve homogeneity. The notion of fractional 
constant is also adressed. Finally, an illustration is given in section [5] with the case of 
natural Lagrangian systems. A complete resolution is given for the harmonic oscillator. 

2. Extrinsic characteristic time 

In this paper, the notation d^/dt" will concern an unspecified fractional derivative, except 
in section 15. 2[ Moreover, the difficulties inherent to fractional derivatives (for example, 
d"" /dt" o d^/dt" 7^ d'^"^ /dt^" in general) won't be addressed here. 

As it has already been said, the fractional aspect of equations occurs with derivatives 
generalized to noninteger orders. If this sustitution is done in one term of an equation, the 
other ones have also to be modified so as to preserve homogeneity. As a first example, we 
consider the damped oscillator equation, 

—x{t) + X—x{t)+uj^x{t) = 0. (1) 
The constants A and a; are homogeneous to T~^. Replacing d/dt by d" /dt" gives 

This equation is non-homogeneous: the three terms have different temporal dimensions 
(T~", T~^^~^°'^ and T~^). A homogeneous version could be 

^x{t) + X-—x{t)+u'-x{t) = 0. (2) 

Things are getting trickier with the diffusion equation, 

d cP' 

—u{x, t) - D-^u{x, t) = 0, with [D] = T-\ (3) 

If the operator d/dt is directly remplaced by 9"/(9t", the only way to preserve the 
temporal homogeneity is to shift D into D". We obtain 

Unfortunately, this equation becomes now nonhomogeneous in space: has no 

spatial dimension, whereas the spatial dimension of D°'d^/dx'^ is The only way 

to preserve the spatial homogeneity is to change d/dx by 0°" /dx'^, which leads to 

^n(a:,t)-D"^n(x,t) = 0. (4) 
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Regarding the numerous studies dealing with that equation (for example, [H[T3]), one 
of the most used in fractional calculus, this case seems too restrictive. According to [I3l 
chap. 16], anomalous transport cannot be described by (Ilj) (because the transport exponent 
remains classical). 

As a last example, we consider the free fall equation, 

^x(t) + (7 = 0. (5) 
The dimension of the gravity acceleration g is [g] = LT~'^. 

If d?/dt^ is changed directly into d'^" /dt^'^, it is impossible to obtain an fractional and 
homogeneous equation. 

A first conclusion could be drawn: a necessary condition for an equation to be physically 
relevant in the fractional case is the possibility of obtaining a homogeneous formulation. 
In this case, the damped oscillator and the diffusion fractional equations fill this condition, 
unlike the free fall equation. 

However, the constants which appear the equations may be built from other constants. 
For example, the gravity acceleration g is in fact defined hy g = Q R^^, where Q is the 
gravity constant, Mt et Rt the mass and the radius of the Earth. So, equation (l5|) can 
provide a fractional and homogeneous equation, 

-^x{t) + g'^R}f- = 0. 

Therefore, the initial equation itself is not suffisant to conclude if the fractional equation 
derived is physically relevant or not. This first conclusion is consequently invalid. 

Moreover, if no restrictions are added, there may be several ways to obtain fractional 
equations, which is unsatisfactory. Indeed, ([I]) can notably turn into 

^<t) + A"^x(t) + u^X'^-Ht) = 0, 

or into 

^<t) + A^"-^^x(t) + u'-x{t) =0, 0<(3<1. 

Both of these equations differ from (l2|). 

Facing this problem of underdetermination, equations have to be constrained so as to 
provide a unique fractional and homogeneous equation. As it was mentioned above, the 
fractional embedding plays this part in the case of Lagrangian systems. It implies, among 
others, that there is a single fractional exponent in fractional equations describing such 
systems. Furthermore, it does not modify the constants of the initial equation. Once the 
Lagrangian is given, the structure is totally determined and the fractional aspect appears 
only through the switch of d/dt into d^ /dt". Consequently, a suitable fractional operator 
should not modify the homogeneity of the initial equation. Two possibilities are conceiv- 
able. The first one is the temporal nondimensionalization of the system: the classical and 
fractional temporal operators become dimensionless, and the homogeneity of the equation 
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is preserved. The second method consists in choosing a fractional operator which dimen- 
sion is the same as the classical one, i.e. T~^. Both of them require the introduction of 
a time constant r. The examples above show that r may not be determined directly by 
the equation (see (El)), or that several constants could suit for r (see ([II)). Consequently, 
the time constant r has to be postulated, and defined a priori, i.e., independently of the 
classical equation. The introduction of such a constant enables to solve the problem of 
homogeneity, while remaining compatible with the fractional embedding. 



3. General embeddings 
We adapt here the embedding notion presented in [1]. 

3.1. Embedding of difFerential operators. Let f = (/i, . . . , /p) and g = {gi, . . . ,gp) be 
two p-uplets of functions M''+^ — > M. Let a,b eR with a < b. We denote by 0(f, g) 
the differential operator defined by 



0(f,g) : C^i[a,b]) — > C^{[a,b]) 



where, for any function u : M'^^^ — > M, any functions Xq, . . . ,Xk ■ [a, b] — > M, 

u{xo{»),...,Xk{»),») ■■ [a,b] — > R 

t I — > u{xo{t),...,Xkit),t), 

and, for any functions / and g, (/ ■ g)it) = f{t) ■ g(t). 

Now we extend this class of operators with an operator which generalizes d/dt. 

Definition 1. Let V : C°°([a,6]) — > C°°([a, 6]) be a differential operator. With the 
previous notations, the V-embedding ofO{f,g), denoted S{0{f, g),V), is defined by 

8{0{lg),V) : C-([a,6]) — > C-([a,6]) 

p 

X ^ 5^(/,-P^o^,) (a;(.),...,P'=x(.),.). 

i=0 

We define the ordinary differential equation associated with C(f, g) by 

O(f,g)(x)=0, xeC^{[aM). (7) 

Definition 2. Let V : C^{\a,b]) — ^ C°°([a,6]) be a differential operator. With the 
previous notations, the V-embedding of equation (jTl) is defined by 

£iOifg),V)ix) = 0, xeC^i[a,b]). 
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3.2. Embedding of Lagrangian systems. We consider a Lagrangian system, with La- 
grangian L{x,v,t), t G [ci,&], and V : C°°([a, 6]) — > C°°{\a,h]) a differential operator. 
The Lagrangian L can naturally lead to a differential operator of the form ([6]) 

0(1,L) : C-([a,6]) 
Now we identify L and 0(1,L). 

The P-embedding of L, £{L, "D), will be denoted L{V) or shortened into L, 

L{V) : C-([a,6]) ^ C-([a,6]) 

X I — ^ 'Da;(«), •). 

In Lagrangian mechanics, the action and its minima play a central role. Let us define 
an action suitable for the embedding formalism. 

Definition 3. Let g : C°°([a, 6]) — > C°°{\a,h]) be a mapping. The action of g, denoted 
A{g), is defined by 

A{g) : C-([a,6]) M 

X I — >• / g{x){t)dt. 



For example, with the identification L = (9(1, L), the action of L is given by 



b 



d 



A{L){x) = J L \^x{t),—x{t),tj dt. 
Concerning the P-embedding of L, the associated action is 

f6 



A{L{V)){x)= I L{x{t),Vx{t),t) dt. 

J a 



The extremum of the associated action to a Lagrangian system provides the equation of 
motion of the system trough the following theorem. 

Theorem 1. The action A{L) is extremal in x if and only if x satisfies the Euler-Lagrange 
equation, given by 

Vt G [a, 6], d^L (^{t), ^a;(t), - j^d^L (^{t), j^x{t),t^ = 0. (8) 

Equation ([8]) is denoted EL{L). For the embedded Lagrangian, a connected result states 
the following. 

Theorem 2. Under some appropriate variations of the action (see \7]), A{L{V)) is ex- 
tremal in X if X satisfies the causal Euler-Lagrange equation, given by 

VtG[a,6], diL{x{t),Vx{t),t) -Vd2L{x{t),Vx{t),t) = 0. (9) 
Equation ([9]) is denoted EL{L(T>)). This leads to the following result. 



6 



PIERRE INIZAN 



Theorem 3. In the causal case (i.e. with the appropriate variations of the action), we 
have the commutative scheme, 



This embedding procedure is qualified as coherent. 

Proof. Equation (El) can be written in the form ((71), with f = (1, 1) and g = {diL, —d2L). 
Its P-embedding is exactly ([9]). 



Let us now explain the two methods presented in the Introduction for the Lagrangian 
sytems. We show that they are equivalent. Then we present what could be a third method, 
also equivalent, which is farther from the philosophy of the embedding theory, but which 
makes fractional constants to appear. 

Let V, t) be a Lagrangian, with t G [a, 6], and r the time constant introduced in the 
previous part. The temporal dimensions of the variables x, v and t are respectively T°, 



4.1. Temporal nondimensionalization. Nondimensionalization is a widespread method 
used in physics, particularly in fluid mechanics, to simplify equations and to exhib relevant 
parameters. Concerning fractional dynamics, the constant r introduced above seems phys- 
ically relevant. In the following part, it is used to nondimensionalize equations according 
to the temporal dimension. 

Definition 4. We note d = a/r and b = b/r. The temporal nondimensionalized La- 
grangian, denoted L^, evolving on the interval [d,b], is defined by 



None of the variables of L„ has temporal dimension, which justify the denomination. We 
choose here V = /dt"". As the new temporal evolution variable u has no dimension, the 
classical derivative according to this new time, d/du, is also dimensionless, such as d'^ jdvP" . 
Therefore, the substitution djdu d^ jdvP" does not break homogeneity anymore. 

For example, if d^ /dt^ is taken as the Caputo derivative, defined by 



we see that if to and t have the dimension of a time, d^'/dt"' = t^Vf is homogeneous to 
T~", whereas if and u are dimensionless, d^/du"" = uo'^u dimensionless. 



EL{S{L,V)) =S{EL{L),V). 



□ 



4. Homogeneous fractional embeddings 



and T. 



VnG[a,6], Ln{x,y,u) = L{x, —,Tu). 




The (i"/(iti"-embedding of L„, L 



'n 



{d°/du'^), is 
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For sake of lisibility, I/„ {d°^ /du") will be shortened into It verifies 

du 



d" 

Vx G C°°([a, 6]), Vm e [a, b], Ln{x){u) = L ( r^-^- — x{u),tu 



The associated action to L„ is 



'a 



A{Ln)ix) = J Ln yx{u),—x{u),uj du 

and leads to the following result. 

Theorem 4. The Euler-Lagrange equation associated with the d" / du"- embedding of the 
nondimensionalized Lagrangian L„ is given, in the causal case, by 

f d" \ d" f d" \ 

\/ue[a,b], ( — j - — a2L„ ( — j = 0, (10) 

with X e C^{\d~b]). 

This is the equation of motion of the Lagrangian system which Lagrangian is and 
which dynamics evolution is fractional, according to the operator d"/du". The solution 
X of ( fTOl l depends on the nondimensionalized evolution variable u. As d"/du" has no 
dimension, equation ( ITOl ) is clearly homogeneous. 

This method eludes the question of the dimension of the fractional operator, and, by 
modifying the Lagrangian variables (but not the Lagrangian itself), it leads to an homo- 
geneous fractional embedding. However, it occults the real dynamics, based on a real (i.e., 
dimensionalized) time. The following method will lift this veil. 

4.2. Homogeneous fractional derivative. The time constant r can be used to nondi- 
mensionalize the temporal operator, but it also enables to build a fractional operator which 
dimension remains T~^: t"~^ d" / dt" . Therefore, this operator preserves the homogeneity 
of the fractional embedding. Let us detail the procedure. 

We consider the same initial Lagrangian L{x,v,t), with t G [a,b]. Contrary to the 
previous method, there is no need to define a new Lagrangian nor a new evolution variable. 
Now the differential operator V becomes T"~^d"/dt" and acts on functions which depend 
on the real time t. 

The T"-^d"/dt" -embedding of L, L {T"-^d"/dt") will be denoted Lh. 

U : C-([a,6]) C-([a,6]) 



X 



The associated action A{Lh) verifies 

Vx G C°°([a,6]), A{U){x) = j\ (^x{t),T"-'^x{t),t^ dt. 



The equation of motion is now given by 
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Theorem 5. The Euler- Lagrange equation associated with the (r° ^ d^ / dt^)- embedding of 
L is given, in the causal case, by 

VtG [a,6], d,L(^x{t),T--'^x{t),t^ -r--'^d,L{x{t),r'^-'^x{t),t^=Q. (11) 

In this method, the fractional Euler-Lagrange equation and its solutions x depend of 
the real time t, and all of the fractional objects (variables, operator, Lagrangian, equation) 
hold the same dimension as their classical predecessors. 



4.3. Equivalence between the two methods. Both of these two methods constitute 
homogeneous fractional embeddings. Fortunately, if a quite natural condition is respected, 
they turn out to be equivalent. 

We suppose that the fractional derivative d°'/dt°' verifies the following condition: 

d" ^ 

Va; e C°°([a,6]), r"— x(t;)|,^, = — x(m)|„^,/, , where x : u ^ x{ut). (12) 

This property is notably verified by the usual fractional derivatives. For example, for 
the Caputo derivative, we have 



aKHu) = J^^ -J {u-v)-'^x\v)dv 

r(l - a) J a 

/ {u--)-^~x\-)-, S = VT 



r(l - Jar ^ ^ ^ 

— / {uT - s)-''{Tx'{s))ds 

r(i -«) Ja 

r(i - a) 

r'^aV^xit). 



{t — s) °'x'{s)ds, t = UT 



The two methods are related through the following result. 

Theorem 6. If the fractional operator d" / dt" respects condition (fT2l ). then the two methods 
are equivalent: x : t ^ x{t) is solution of (fTTI) if and only if x : u t-^ x{ut) is solution of 

Proof. Partial derivatives of L and L„ verify 

diLn{x,y,u) = diL (^x,^,UT^ , (13) 
d2Ln{x,y,u) = ^d2L (^x,^,uT^ . (14) 
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From (111, (111 and (fll, we infer 



^ \ ( d"' 

diLn[x{u),—x{u),u\ = diL lx{t),T''-^—x{t),t 

^ d°' _ \ If d°' 

d2Ln[x{u),—x{u),uj = -(92^ ( a;(t),r""^— x(t),t 



which proves that (jlOl) and (jlTi ) are similar. 

□ 

Therefore, if condition (fT2l ) is verified, nondimensionalized solution u E [a/r,b/T] x{u) 
of (fTOll can be redimensionalized into t e [a, 6] ^ xit/r), solution of the dimensionalized 
and equivalent equation (fTT]l . 

4.4. Fractional constants. In various litterature dealing with fractional equations using 
the operator d°'/dt°' (homogeneous to T~"), homogeneity is sometimes preserved by in- 
troducing fractional constants (see [9] for example). For a wide class of Lagrangians, the 
(i"/(it"-embedding of a modified Lagrangian leads to a homogeneous fractional equation 
involving fractional constants. 

We consider the same Lagrangian system as above. The temporal dimension of L is here 
denoted T"" (uq = —2 in most of the cases). We suppose that L is equal to its Laurent 
series in the variable v, 

L{x,v,t) = ^aifi{x,t)v\ (15) 

where [a,] = T*"'"'^°. The functions fi have hence no temporal dimension. 

Using the fractional operator T°'~^d°'/dt"', we obtain for any function x G C°°{[a,b]) 

= ga,/,(x,t)r^("-^) [^x{t)y 
= ^^ifi{x,t)[-^xi't)] , with 



The dimension of di is [di] = T*""^"". With this fractional dimension, those coefficients 
can be considered as fractional constants. 

We can form a new Lagrangian with those constants. 

Definition 5. For a Lagrangian L{x,v,t) = '^i^iCLi fi{x,t) v\ the fractional Lagrangian, 
denoted Lf, is defined by 

Lf{x,v,t) = J^ai/i(x,t)' 
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The (i"/(it"-embedding of Lj, denoted Lf, verifies 

Lf : C-([a,6]) C^i[a,b]) 



Hence we have 



which leads to 



Theorem 7. The Euler-Lagrange equation associated with the d'^ /df^- embedding of Lf is 
given, in the causal case, by 

.c^ / V / 

withx e C°°{[a,b]). 

With relation (fT6l) . we infer the following result. 

Theorem 8. // condition (fTSl ) «5 verified, the {t"~^ d^ / dt")- embedding of Lf and the 
d" / df^- embedding of Lf are equivalent, 

In the causal case, we have besides ffTSll = ffTOD. 



One verifies that equation f flTl ) is homogeneous: inhomogenity of the operator d'^/dt"' 
is counterbalanced by the fractional constants a,. However, the homogeneity problem is 
here only shifted: Lf {x{t),d/dtx(t),t) and EL{Lf) are now inhomogeneous. Furthermore, 
the physical relevance of Lf is discutable, and this Lagrangian may appear as an ad hoc 
construction to obtain homogeneous fractional equations. In this sense, this third method 
does not fit with an underlying idea of the embedding theory. For any Lagrangian system, 
its Lagrangian is independant of the dynamics: it is only the choice of the evolution 
operator, i.e. the temporal derivative, which conditions the dynamics of the system. 

Finally, it seems that the second method is the most relevant homogeneous fractional 
embedding. 

5. Natural Lagrangian systems 



5.1. General case. An illustration of those methods is given here with the most wide- 
spread type of Lagrangians: the natural ones. 
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Definition 6. We consider an autonomous Langrangian system with Lagrangian L{x,v), 
evolving on a temporal interval [a, b]. L is a natural Lagrangian if it can be written as 

L{x, v) = ]-mv'^ - V{x), (20) 

where m is the mass of the system, |mf ^ is the kinetic energy, and V{x) the potential 
energy. The temporal dimension of L is T~^. 

We introduce a time constant r, and we apply now the three embeddings to those 
systems. 

5.1.1. Temporal nondimensionalization. We denote a — a/r and b — bjr. 
The new Lagrangian is 



Ln{x, y)^L (x, ^) = ^my'^ - V{x), 
and its embedding verifies 



1 / (i" ^ ^ 



yxeC^{[d,b]),yue[a,b], Lr,(x){u)^—m\^—x{u)j -V(x(u)). 
The associated action is 

A{Ln)ix) ^ J Ln (^x{u), ^^(^)^ du 

and tlie Euler-Lagrange equation, in the causal case, turns into 

5.1.2. Homogeneous fractional derivative. The initial Lagrangian remains inchanged. Its 
embedding is now 

VxeC~([a,6),Vte[a,6], Ln{x){t) = ^-rnr^^"-'^ (^^x{t)^ -V{x{t)), 
the action is given by 

A{Ln){x) = j\(^x{t),T"-'^x{t),t^ dt 

'r^^'^'' (^<t)y -V{x{t))dt, 



and the causal Euler-Lagrange equation is 



J2a 

Vte[a,6], mT'i--^^^x{t) + V'{x{t)) = 0. (21) 
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5.1.3. With fractional constants. Let xq be a position such as V{xo) ^ 0. The potential 
energy can hence be written as V{x) = aofo{x,t), with oq = V{xo), which temporal 
dimension is T~^, and fo{x,t) = V{x)/V{xo), which is dimensionless. 
The Lagrangian (i20l ) can therefore be expressed in the form (fTSl l. with 

• aj = for i < 0, 

V(x) 

• ao = V{xo), fo{x,t) = 

• ai = 0, ^ 

• a2 = m, f2ix,t) = -, 

• ttj = for j > 2. 

The fractional constants are 

• flg = ao = V{xo), 

and lead to the new fractional Lagrangian 



Lf{x,v) = ^mr^(""^)i;2 - V{x). 



The associated action is defined by 

AiLf)ix) = flrnr'^"-'^ (^^^^0' ^ V{x{t))dt. 
Its extremum provides the following Euler-Lagrange equation: 

Vt G [a, 6], mT'(--''>^x{t) + V (x(t)) = 0, 
which is similar to (l2Tl) . 

Now, let us solve those equations in the case of the harmonic oscillator. 

5.2. Harmonic oscillator. This system has already been studied in the fractional case 
(see [mini [12]). 

The harmonic oscillator is a natural Lagrangian system which Lagrangian L is given by: 

L{x,v) = -mv — -kx , (22) 

with [x] = L, [v] = LT-\ [m] = M et [k] = MT-^. 

As above, we introduce a time constant r (which is a priori not linked with ^Jrnjk). 
We arbitrarily choose here the Caputo derivative for the fractional operator /dt"'. 

5.2.1. Temporal nondimensionalization. The nondimensionalized Lagrangian L„ is: 

1 1 
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Its embedding, Ln{x){u) = ^m{d°'/du°' x{u))'^ — ^kx{u)'^ leads to the following Euler- 
Lagrange equation: 

E \a,b], -7; f n x(u) + kxiu) — 0. 



With the usual notation uj — ^Jk/m, we obtain 

^ + (wt)^x(m) = 0. 

The (nondimensionalized) solution, denoted is given by 

e [a, 6], = Xn{d)E2a {-{ujrf{u - df") , 

where, for A > 0, is the Mittag-Leffler function defined on C by 



,.„r(Afc + i)- 

The dimensionalized solution, denoted xa, is hence 

Vt G [a, b], Xd{t) = Xn{^) = Xd{a)E2^ {-uj^r^^^'^^t - a)^") . 
One verifies that the argument of the Mittag-Leffler function is dimensionless. 

5.2.2. Homogeneous fractional derivative. The embedded Lagrangian Lh, given by 

2 



provides the equation 



4 = imr^("-) (l^x(t)) -kx{t)\ 



J2a 

e [a, 6] , mr'^'^-^) ^x(t) + kx(t)^0. 



(23) 



Introducing a;, we obtain 

yte[a,b], ^a;(t)+a;V(^-")a;(t) = 0, 
which admits the solution Xh defined by 

Xhit) = xu{a)E2a (-cuV(^-")(t - af-) . 
One verifies that Xn — Xh. 
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5.2.3. With fractional constants. The Lagrangian (l22l ) can be written as (fTSll . with 

• ao = k, fo{x,t) = 

• ai = 0, 

• 02 = m, /2(x,t) = |, 

• = for i > 3. 

The fractional constants are given by 

The associated fractional Lagrangian, -Z^/, is 

We obtain the Euler-Lagrange equation, 

VtG[a,6], mr2("-i)^x(t) + = 0, 

which is similar to equation fl23l) . 

5.2.4. Comparison with the inhomogeneous fractional embedding. If we directly apply the 
fractional embedding presented in [4] to the Lagrangian (l22l l. we obtain 

^x(t)+^2a;(t) = 0. (24) 
The temporal dimension of d"^" / dt^'^ x{t) is T"^", whereas for to'^x(t), it is T"^. 



The solution Xj of (|2lD is 

The argument of the Mittag-Lefiler function is whereas it should be dimension- 



less. 



6. Conclusion 

With the first two embeddings presented above (temporal nondimensionalization and 
homogeneous fractional derivative), it becomes possible to preserve homogeneity of em- 
bedded equations. Both of them lead to the same fractional equation. Therefore, from 
any Lagrangian system, it is possible to construct a fractional equation which respects two 
fundamental physical principles: homogeneity and least action principle. 

To do so, a time constant r has to be introduced. Similarly to the fractional exponent a, 
this parameter is not defined by the equation itself. Contrary to the nondimensionalization 
of classical equations, fractional solutions depend on r (see the example of the harmonic 
oscillator). This time constant is hence more than a calculus intermediary: it conditions the 
dynamcis of the system. Its physical relevance seems as much important as the constants 
of the initial equation. To sum up, a fractional equation should be characterized by two 
parameters: the fractional exponent a and the time constant r. 
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So as to understand the physical meaning of this characteristic time, it may be useful 
to set the context in which such equations arise. In [5l [I3], the fractional aspect of the 
dynamics appears only on long time scales. At a microscopic level, the dynamics is classical 
and evolves according to a microscopic time, but when one zooms out at a macroscopic 
level, ruled by a new macroscopic time, the fractional aspect may emerge. It might be 
through the relation between these two time scales that the time constant r may appear. 
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